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numerical analysis on the performances of a SUPPS identical to the prototype in Manzanares, Spain 
which is exposed to the external crosswind with different velocities is carried out in this paper. 
A geometrical model including the SUPPS and its outside ambience is built and the mathematical 


Keywords: models to describe the fluid flow, heat transfer and output power of the whole system are further 
Solar updraft power plant system (SUPPS) developed. The pressure, temperature and velocity distribution of the air in the ambience and SUPPS 
A crosswind together with the output power of the SUPPS are analyzed. The numerical simulation results reveal that 
Collector ambient crosswind has influence on the performance of the SUPPS in two ways. On one hand, when the 


ambient crosswind is comparably weak, it will deteriorate the flow field and reduce the output power 
of the SUPPS. On the other hand, it may even increase the mass flow rate and output power if the 
crosswind is strong enough. 
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Nomenclature 


Cz1,Ce2,Cz3 Constants for turbulent model 


Cp specific heat at constant pressure (J/(kg K)) 

g acceleration of gravity (m/s?) 

G solar radiation intensity (W/m?) 

Gk turbulence kinetic energy generation due to the mean 
velocity gradients (J) 

Gp turbulence kinetic energy generation due to 
turbulence (J) 

L collector height (m) 
heat flux through the ground underneath the collec- 
tor (W/m?) 

Ra Reynolds number (dimensionless) 

T temperature (K) 

t time—for unsteady items (s) 

u velocity in y-direction (m/s) 


U200m ambient crosswind velocity at the height of 
200 m (m/s) 

w velocity in z-direction (m/s) 

Xyz Cartesian space coordinates 


a thermal diffusivity (m?/s) 

v kinetic viscosity (m?/s) 

B volume coefficient of expansion (1/K) 
p density (kg/m?) 

T shear stress caused by viscosity (N/m?) 
K karman constant 


Subscripts 
s surface value 
ij any direction of x, y and z 


They copy the daily solar thermal air motion in the atmosphere to 
produce electric energy free of CO emissions, and are predicted 
to be an efficient way of mitigating the unprecedented pressure to 
reduce CO2 discharge that many countries in the world are facing 
today. Compared with more conventional solar energy applica- 
tions, the SUPPS has gained its ascendancy by achieving several 
goals: easy to procure building materials, less contaminants 
generation throughout its operating process and a longer operat- 
ing life span. A SUPPS often consists of four key parts: a collector, 
a chimney, a turbine and an energy storage layer. The collector, 
whose canopy is made of transparent or half-transparent materi- 
als, such as plastic and glass, is large enough to collect solar 
energy because of the one-way-screen characteristic of the 
canopy material. During the day the visible and UV wavelengths 
of the solar radiation, where most energy of the sunlight spec- 
trum is concentrated, pass through the cover and warm up the air, 
meanwhile the infrared wavelengths warm up an energy storage 
layer. The increased air temperature results in a decrease of its 
density, while in the meantime a strong airflow buoyant force is 
produced and the chimney stack leads to natural convection in 
the SUPP. The air momentum will drive a turbine at the foot of 
solar chimney stack whereby kinetic energy will be transformed 
into electric power. During the night, the heated energy storage 
layer being made of soil, stone or water in tubes, transfers 
thermal energy to the air inside the collector, allowing night 
operation. 

The first 50-kW SUPPS prototype built that lead to the 
validation of the solar chimney concept was originally erected 
in Manzanares, Spain by Schlaich [1], a professor at University of 
Stuttgart, in the early 1980s as a result of a joint venture between 
the German government and a Spanish utility. The chimney was 
194.6 m high and 10 m in radius, with the collector being 122 m 
in radius and 2-6 m in height from inlet to center. This prototype 
power station worked successfully for more than 7 years. 

Since then, relevant studies on the SUPPS have never ceased 
due to its prosperous future and some significant breakthroughs 
in theoretical, numerical analysis and prototype experiments 
which have been since carried out. All around the world several 
teams are examining the possibilities of building SUPPS for 
instance Enviromission in Australia and in Arizona (USA) [2,3]. 
Larbi et al. [4] studied the possible performance of a SUPPS in the 
south western region of Algeria, and Dai et al. [5] in the north 
western regions of China. Zhou et al. [6] did the same for the 
Qinghai-Tibet region of China, Mostafa et al. [7] estimated the 
performance of a solar chimney under Egyptian weather 


conditions, Sangi [8] for Iran and Ketlogetswe et al. [9] studied 
the case of Bostwana, Hamdan [10] the Arabian Gulf region, 
Nizetic [11] looked at the feasibility of implementing SUPPS in 
the Mediterranean region, as well as Cervone et al. [12] and Bilgen 
and Rheault [13] developed a mathematical model to evaluate the 
performance of SUPPS at high latitudes. 

Based on the 50-kW prototype in Manzanares, Spain, Haaf 
et al. [14,15] made primary investigations into the energy 
balance, design criteria, and cost analysis in the SUPPS. Following 
that were the same author’s experimental reports on the operat- 
ing condition of the SUPPS in Spain. In order to analyze the 
influence of miscellaneous parameters, such as environmental 
conditions and geometrical dimensions, on the temperature and 
velocity of air and output power of the solar chimney, the 
research group led by Sherif [16-21] conducted comprehensive 
mathematical models to evaluate the fluid flow, heat transfer, and 
output power performances of various scales of SUPPS and 
developed three types of experimental prototype in Florida, with 
the chimney shape, collector construction and energy storage 
layer performance being taken into consideration. 

By calculating the performance and efficiency of a SUPPS with 
chimney friction, turbine, and kinetic energy losses being con- 
sidered, Gannon and von Backström [22,23] brought forth an air 
standard cycle analysis of the solar chimney power plant, accom- 
panied by more thorough analysis of SUPPS with turbines being 
employed by Kroger and Buys [24], Gannon and von Backstrom 
[25,26] and Ming et al. [27]. Bernardes et al. [28-30] developed a 
comprehensive mathematical model to analyze large scale SUPPS 
with a double and single collector canopy with an energy storage 
layer and turbine performance being considered, comparing 
simulation predictions to experimental results from the prototype 
plant at Manzanares, and evaluated the operational control 
strategies applicable to SUPPS. Schlaich et al. [31] presented the 
basic theory, practical experience, and economy of SUPPS to give a 
guide for the design of 200-MW commercial SUPPS. Pretorius and 
Kroger [32,33] developed a comprehensive mathematical model 
to numerically simulate the SUPPS and analyzed the impact of 
different calculating methods on output power. 

Ming et al. [34] developed a comprehensive model to evaluate 
the performance of a SUPPS in which the effects of various 
parameters on the relative static pressure, driving force, power 
output and efficiency were further investigated. Zhou et al. [35], 
Kasaeian et al. [36] and Ferreira et al. [37] conducted various 
experimental analyses on mini-scale SUPPS. Koonsrisuk and 
Chitsomboon [38,39] and later Sangi et al. [40] performed 
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detailed theoretical and numerical simulations of SUPPS. Maia 
et al. [41] theoretically evaluated the influence of some 
parameters on the behavior of the airflow in a solar chimney 
indicating that the height and diameter of the chimney are the 
most important physical variables. Ming et al. [42] and Zheng 
et al. [43] further studied the power generation and efficiency of 
the solar chimney power plant systems coupled with turbine and 
discussed effects of different numbers of blades in the Spanish 
prototype. Also a MW-graded SUPP model, followed by the 
numerical analysis the influence of energy storage layer on the 
fluid flow, heat transfer performances of SUPPS by Ming et al. [44] 
and Xu et al. [45]. As well as this, Ming et al. [46] conducted a 
numerical analysis on the selection of chimney shape and 
chimney ratio for a 10 MW SUPPS with the aim of achieving the 
maximum output power with minimum cost. The results indicated 
that the cylindrical chimney would be the best choice among the 
three basic configurations (divergent, conical and cylindrical chim- 
ney), whose optimum H/D value ranges from 6 to 8. 

In the latest four years, quite a few research work on SUPPS 
has been published. Fluri and von Backström [47,48] compared 
the performance of different turbogenerator layouts, single rotor 
and counter rotating turbines, both with or without inlet guide 
vanes, using analytical models and optimization techniques, and 
discussed the important design parameters. Koonsrisuk and 
Chitsomboon [49-51] predicted the performances of large scale 
SUPPSes using the dimensional analysis together with engineer- 
ing intuition to combine eight primitive variables into only one 
dimensionless variable that establishes a dynamic similarity 
between a prototype and its scaled models. Zhou et al. [52] 
analyzed the influence of chimney height on the performance of 
SUPPS subjected to standard lapse rate of atmospheric tempera- 
ture. Zhou et al. [53] and Fluri et al. [54] conducted detailed 
economic analysis on SUPPSes. Petela [55] conducted a simplified 
interpretative mathematical model of the SUPPS which could be 
used to demonstrate feasibility of application of exergy for 
analysis of SUPPS and for proposing the methodology of the full 
thermodynamic analysis including exergy. This study also pre- 
sented the application of the concepts of exergy and gravity input 
for the modified exergetic interpretation of processes. Reports on 
different applications of conventional SUPPS combined with other 
systems and sloped SUPPS can be found in [56-61]. 

However, a close look at scientific publications shows that 
researchers focused more on solar effects than in crosswind 
influence outside the solar chimney. Yet it is generally accepted 
that, the influence of ambient crosswind on the performance of 
the SUPPS is self-evident, even as significant as the influence of 
the solar radiation. Up to now, only a few preliminary studies 
have been carried out. 

Many authors like Niemann and Hoffer [62] were mainly 
interested by structural or architectural aspects, concerning 
mechanics and concrete resistance to wind, vibration or earth- 
quakes. A model with crosswind was studied by Niemann et al. 
[63], but was only limited to structural reliability of SUPPS 
system. Rousseau [64] showed that structural integrity of solar 
chimneys might be compromised by the occurrence of resonance. 
The wind gust spectrum peaks near the solar chimney’s funda- 
mental resonance frequency which poses a reliability threat, not 
only to the solar chimney, but also to all high-rise, slender 
structures. van Zijl and Alberti [65] presented the results of series 
of physical experiments in wind tunnels establishing external and 
internal pressure coefficient distributions and overall drag coeffi- 
cients for rigid smooth cylinders, and demonstrated the stabiliz- 
ing role of rigid cylinders with vertical ribs. Harte and van Zijl [66] 
studied the static wind profile converted to pressures acting on 
the chimney along its height, as well as along the circumference. 
Lupi [67] and Borri et al. [68] performed innovative modeling of 


dynamic wind action on SUPPS, and studied structural optimiza- 
tion of solar towers to minimize wind induced effects but again 
limited to structural dynamics aspects. Harte et al. [69] showed 
that natural draft cooling towers and chimneys of SUPPS have 
many structural properties in common: they are shell structures 
made of reinforced concrete, they transport by their internal 
updraft warm air into the atmosphere, and because of their 
height, gale actions play the most important role in the design 
and show how far structural design problems of these structures 
are common. Kraetzig et al. [70] explained the design of high 
efficient tower shells for SUPPS and CT including their critical 
response characteristics and demonstrated their close structural 
mechanical relationship to each other. Lv et al. [71] studied SUPPS 
structure vibrations and frequencies. 

Already in the 1970s large dry cooling towers reaching up to 
300 m, were designed for power stations in arid zones. However, 
at this height they was a long way off the height needed for 
professional competitive and operational solar chimneys, which 
start at approximately 700 m. 

Although the operation of solar chimneys and cooling towers 
(CT) are quite different, the problems arising from wind gusts and 
reduction of air-intake-flow-rate under crosswind conditions 
which decrease the operation efficiency are well known because 
of available real scale results obtained from the large number of 
CT built over the last 40 years. For SUPPS there is little data 
available, as the tallest prototype 195 m high was built in 1982 in 
Manzanares (Spain), and since then one smaller one has been 
built. The promoters of the Arizona solar tower [72], and of the 
Wuhai Inner Mongolia SUPPS [73] benefit from little and scarce 
real scale experience and data to strengthen their preliminary 
studies. The efficiency of SUPPS depends mainly on the size of the 
collector area and on the height of the chimney, both reasons for 
their enormous dimensions: collector diameters up to 7 km and 
chimney heights up to 1500 m are on pre-design. 

As the concrete cement has to prove its realistic performance 
by enduring resistance under environmental and operational 
conditions, constructing cooling tower shells from 100 to 300 m 
high is a big challenge and the wind load has been well studied. 
The challenge is enhanced tenfold for solar updraft towers from 
700 to 1500m high: high compression stresses under dead 
weight and wind action, suction or forced wind vibrations in 
the upper chimney part. The construction of all these very tall 
chimneys will make high demands on construction techniques 
and concrete technologies. 

But even if CT have longer seniority, till recently the phenom- 
ena’s responsible for wind decreased performance was not well 
identified and has grabbed the attention of researchers for many 
years. The way the CT operates under crosswind conditions is still 
at the forefront of energy research, whereas in the area of SUPPS it 
is still in its infancy. 

For CT, earlier studies made in 1988 by Radosavljevic and 
Spalding [74], simulated the fluid flow and temperature distribu- 
tion in a wet cooling tower affected by cross-wind. The wind 
effect on a natural draught dry cooling tower was studied in 1993 
and 1995 by Du Preez and Kroeger [75,76], using a simple 
turbulence model based on the eddy-viscosity to study the flows 
inside and on the periphery of the tower. In order to minimize the 
crosswind effect on the thermal performance of towers, they 
introduced the windbreak walls. Wei et al. [77] used in 1995 
complete sample models and wind tunnel testing to demonstrate 
the adverse effects of wind on the functioning of dry cooling 
towers. 

Problems encountered in the operation of cooling towers in 
the west of Canada conducted Derksen et al. [78], to investigate in 
1996 the effects of wind on the air intake flow rate. Wind tunnel 
tests on a 1/25 scale model allowed the study of the external flow 
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patterns, pressure characteristics, and air intake flow rate. They 
discovered that flow imbalance was the cause of ice formation on 
the inlet side, as the wind-facing side of the tower presented a 
flow rate almost increased by 45% meanwhile the opposite rear 
side flow rate decreased by 19%. 

The numerical simulation performed by Su et al. [79] 
explained the main reason for decline of the performance of 
dry-cooling tower under crosswind and provided some founda- 
tion to improve its thermal performance under cross wind 
conditions. 

Using the commercial software Fluent, in 2004, Al-Waked and 
Behnia [80], did a 3D modeling of a tower and simulated the 
internal and external flows introducing the wind speed profile as 
a new important parameter, to correctly simulate the flow around 
the tower. They also suggested the windbreak walls as a solution 
for reducing the adverse effects of wind. In 2006, Zhai and Fu [81] 
performed a numerical and experimental study of the effects of 
two windbreak sidewalls on the efficiency of a dry cooling tower. 
In 2010, Al-Waked [82] investigated the effects of crosswinds on 
the thermal performance of natural draft wet cooling towers. 

Goodarzi [83] demonstrated that a natural draft dry cooling 
tower is significantly affected under crosswind conditions and 
performance might decrease up to 75% in the range of moderate 
to high wind-velocity conditions. He proposed a new exit config- 
uration for tower stack that could reduce the throttling effect of 
deflected plume and increase cooling efficiency. 

Very little work has been published about crosswind effects on 
the performance of SUPPS: Serag-Eldin [84] addressed the degra- 
dation of SUPPS performance taking into consideration external 
ambient crosswind, and introduced the concept of controllable 
flaps aimed to reduce the proportion of hot air gone with the 
wind by-passing the chimney stack, but did not continue his 
analysis further on effects of different ambient wind profiles; 


velocity inlet symmetry surface 


ground wail 
wall of canopy 


wall with a heat flux 


pressure outlet 


wall of chimney 


Pretorius and Kröger [85] analyzed the influence of ambient 
winds at 2 m/s by regarding its effects on the annual power 
output of the solar chimney power plant. The comparison of two 
models which were simulated gave with quiet ambient condi- 
tions an annual power output of 373.2 GWh and with windy 
ambient conditions 336 GWh (a drop of approximately 11%). The 
windy conditions resulted in an increased convective heat trans- 
fer coefficient, facilitating a greater heat flux from the collector 
roof to the environment and ultimately lower power output. 

In this paper, we numerically analyze the influence of various 
magnitudes of ambient crosswind on the pressure, velocity and 
temperature distributions, heat transfer and output power per- 
formance of SUPPS. 


2. Model description 
2.1. Geometric model 


In this paper, a simplified model of the SUPPS Manzanares 
prototype [1] is adopted for the numerical simulation. As shown 
in Fig. 1, the model has a 200-m-high and 5-m-radial chimney 
and a collector which covers the earth surface in a round shape: 
120m in radius and 2m in height. In order to simulate the 
performance of the SUPPS exposed to such vast space, presum- 
ably, we place the model in the center of an actually nonexistent 
cubical box with its lengths in x, y, z directions of 400 m, 400 m 
and 300 m respectively. A box that has every surface set with 
different boundary conditions which will be displayed in the 
ensuing introduction. In this figure, x-axis is aligned in the 
velocity direction of the ambient crosswind, z-axis in the straight 
up direction. Assuming the symmetric property to be perpendi- 
cular to the y-axis direction, only half of the whole system can be 


pressure oulet 


Fig. 1. 3-D geometrical model including SUPPS and ambience and their boundary conditions. 
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taken into consideration as shown in Fig. 1, thus the whole 
geometric dimensions of cubical box are 400m, 200m and 
300m in x, y, z directions, respectively. Furthermore, in this 
model, the influence of energy storage layer was not considered 
and the geometrical model was not included in the computational 
model as the present PC could not tolerate the number of meshes 
of the model including the ambience and the SUPPS with the four 
parts (collector, turbine, chimney, energy storage layer). 

By comparison, the simulation with a much longer X value, say 
600 m, downstream of the chimney is also performed and the 
influence of this length downstream the chimney on the accuracy 
of numerical simulation results has been presented in Figs. 3 and 4. 
The main aim of this research is to analyze the influence of ambient 
crosswind on the performance of SUPPS. As for the influence of the 
outflow from the chimney outlet on the local ambient climate, Zhou 
et al. [86] has presented some relevant numerical results, and may 
be one of our further study points will also focus on this. 


2.2. Mathematical model 


For the SUPPS, we consider the fluid flow inside is natural 
convection induced by solar radiation heating the ground wall. 
A criterion number that can measure the intensity of the buoy- 
ancy-induced flow is the Rayleigh number being defined as 
follows: 


_ gBATL? 
~ W 


Ra (1) 
where AT is the maximum temperature increase within the 
SUPPS. a, f, and L are the thermal diffusivity, thermal expansion 
coefficient and the collector height, respectively. After prelimin- 
ary calculation, we found that Ra is higher than 10!° which is the 
upper limit transition value of natural convection from laminar to 
turbulent flow, therefore the turbulent mathematical model need 
to be selected to describe fluid flow within the SUPPS. Also, the 
density variation in the whole computational model including the 
SUPPS and ambience is so small that it can even be neglected in 
Mass and Energy equations, making it necessary to resort to 
Boussinesq approximation without causing large errors [8]. As a 
result, Mass equation, Navier-Stokes equation, Energy equation 
and standard x-¢ equations can be written as follows: 
Continuity equation: 


Op apu) _ 
at OX; ais (2) 
Navier-Stokes equation: 
apui) _ Wpujuy) p Ny 

at ay O ax, O @) 
Energy equation: 
ApCpT) , 8(pCpujT) ð (,eT\ Ou; , ôp op 

t ta ax \Aox,) Ha, Tate) O 
Equation for the turbulent kinetic energy k: 
ô ð ð ok 
= (pk)+ ax; (pkui) = x; (cutter x) + Gy, + Gp—pE—Yu + Sp (5) 


Equation for the energy dissipation: 


ô a a ðe E€ e 
aP ax; P5 ax (teer a) F Cie g (Gk + C3cGp)—Coop q Re +Se 


(6) 


where G, denotes the generation of turbulence kinetic energy 
because of the mean velocity gradients and can be defined as 
Gk = —pu;u; (Ou; /3x;). Or, Ck and c; represent the turbulent Prandtl 
numbers for T, k, and e respectively: or=0.9, ok=1.0, o,=1.3 and 


cı and c2 are two constants for turbulent model: c,=1.44, 
C2=1.92. u, = (Cupk? /e) and c,,=0.09. 

The reason why we select the standard «-¢ equation is that the 
authors simplified the turbine installed at the base of the chimney 
as a reverse fan as shown in Section 2.3. The complex geometrical 
construction of the three dimensional turbine model within the 
chimney base will result in complex strong vortexes and turbu- 
lent flow which could be accurately simulated using RNG or 
Realizable x-¢ model [87]. However, in this work, the turbine 
model is treated as a reverse fan: two dimensional boundary 
condition which will not result in complex strong vortexes and 
turbulent flow, so the standard x-e equation is accurate enough 
to describe the present problem. Detailed numerical simulation 
coupled with three dimensional turbine will be performed in the 
future study. 


2.3. Boundary conditions 


When taken into account the ambient crosswind on the SUPPS, 
boundary conditions for both the SUPPS which have been shown 
in the previous publications and the ambient should be carefully 
given. Fig. 1 also represents the boundary conditions set in this 
case as well as the coordinate applied to the model. Detailed 
descriptions of the boundary conditions are shown as follows. 


(1) Inlet boundary(surface at x=0) 

The cases in this paper are based on the assumption that the 
ambient crosswind is fully developed and the temperature 
constantly at 293 K before flowing into the internal space of 
this model. According to the logarithmic law of the wind 
speed profile in atmospheric boundary layer which was 
proposed by Prandtl in 1932 the ambient wind inlet velocity 
can be fitted as the equations below [88]: 


v=w=0 (7) 


u=1/K(ts/p)'/*In(z/Zo) (8) 


in which t, stands for the ground surface shear stress and Zo 
for aerodynamic roughness length of the ground. Specific data 
of Zo for different terrains is available in charts such as the one 
that was collected by Cermak [88]; and in this case x and Zp 
are chosen as 0.4 and 0.01 m, respectively, since a flat desert 
terrain type is pre-set. After that, ts can be calculated from a 
given value of wind speed u at a known height, a height in this 
case is selected that of the tip surface of the SUPPS. 
Outlet boundary (surfaces at x=400 and z=300) 
After influences of the ambience crosswind and the chimney 
outlet being taken into account, there are two outlet bound- 
aries in the model as shown in Fig. 1: the top surface of the 
model for the exit of the wind from the solar chimney and the 
surface of the box parallel to that of the inlet ambient 
crosswind. Pressure outlet boundary condition is applied to 
each of these and the simplification that reversed flow is 
normal to the boundary surface. 
(3) Ground boundary(surface at z=0) 
Both the ground underneath the canopy of the collector and 
the region outside are contained in the ground boundary 
condition. The exposed ground outside the collector is 
assumed to be isothermal and is set to be a wall with a 
temperature of 318 K. This coarse assumption may be more or 
less influence the accuracy of the simulation results compared 
to the practical working conditions outside the SUPPS but will 
not significantly influence the main object of this paper. 
Consequently, functions for boundary layer are employed 


~ 
N 
Vw 
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here in order to simulate the near-surface flow: 


Us = Vs = Ws = 0 (9) 
Ts = plK/In(z/Zo)}° (U? +v?) (10) 
k= ts/(pC;/*) (11) 
&=(ts/p)/? /Ky (12) 


where the subscript “s” represents “surface value”, while 

other values without such a subscript are evaluated at nodes 

that are near the ground surface. According to Eqs. (11) and 

(12) the values of k and e are fixed respectively in these zones. 

The ground surface covered by the collector canopy is 

assumed to generate heat fluxes of different values according 

to various solar radiation intensities. And other parameters 

set for the wall boundary are quite similar to those of the 

ground outside the collector. Such simplifications ascribe to 

the assumptions below: 

a) the vibration of wind effects is neglected and thus the time 
term is considered steady; 

b) the solar radiation is uniform vertical incident rays; 

c) the ground layer is homogeneous and isotropic; 

d) local heat equilibrium has been achieved between the 
ground and the air that bypasses it; 

e) radiation heat transfer among the walls of SUPPS in the 
model is negligible. 


All the assumptions introduced are aimed at avoiding analyz- 
ing more than one factor simultaneously without deviating 
too much from the real conditions. 

Side wall (surface at y=200) 

Due to the long distance of side wall from the solar chimney 
inlet and outlet, the influence of the side wall is relatively 
small compared with others listed above. There is hardly any 
amount of kinetic or thermal turbulence here and thus heat 
and mass transfer between the geometric model and the 
outside ambience is weak on the side wall. Accordingly, the 
parameters of the properties of side wall are chosen to be 
default during numerical simulation. 

Symmetry surface(surface at y=0) 

As mentioned above, because the model in this case is 
symmetric in y-direction, and so is the inlet crosswind, 
symmetry boundary condition is set in this case in order to 
alleviate the computing process. As is shown in Fig. 1, the 
computing length in y-direction is shortened to just its half 
due to the symmetry scheme, thus reducing the computing 
grids in a large scale. The only necessary procedure to get a 
whole-zone field is to mirror the simulated velocity and 
temperature field to the other side of the symmetry surface. 


Fig. 2. Grid distribution of the geometric model 
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(6) Turbine coupling 

3-D Numerical simulation of the SUPPS couple with turbine 
conducted by Ming et al. [42] indicated that it is a little 
difficult to simulate the turbine region and much more 
meshes are needed to accurately describe flow, heat transfer 
and output power performances of the system. For the 
numerical model in this paper it is impossible to realize the 
simulation procedure simultaneously including regions of the 
SUPPS, the ambience and the 3-D turbine due to the limita- 
tion of grids number. Fortunately, however, the research work 
conducted by Pastohr et al. [89] indicated that it is also an 
efficient way to realize the object by simplifying the 3-D 
turbine to be a 2-D reversed fan with pressure drop across it 
being pre-set. This method was verified by Xu et al. [45] and 
Ming et al. [46] and was proven to be effective to alleviate the 
mesh pressure by 3-D turbine region without significantly 
total performance of SUPPS. Thereby, the turbine is regarded 
as a reversed fan with pressure drop across it being pre- 
set although 3-D model for the SUPPS and the ambience is 
selected. To simulate different output power of the turbine, 
we assign the pressure drop a group of values ranging from 
0 to 200 Pa at an interval of 20 Pa. We can thereby calculate 
the output power of this model according to the equation 
below: 


W.=n,ApV (13) 


where We stands for the electric output power of the turbine, 
N, represents the overall energy conversion efficiency from 
thermal to electricity, which is preset as 0.72. This number is 
the multiply of efficiency from thermal to turbine shaft 
output power and that from turbine shaft output power to 
electricity, the former being about 0.8 by Schlaich et al. [31] 
and the latter being 0.9 which is available easily. 


2.4. Meshing skills 


In general, for the same meshing zone, hexahedral (HEX) 
meshing method is more economical and can reduce false diffu- 
sion more efficiently than tetrahedral one. As a result, HEX grids 
were applied in the model and the mesh generation procedure of 
the whole geometric model was executed using the commercial 
software package Gambit 2.3.16. Fig. 2 reveals the grid distribu- 
tion of the geometric model in different angle of view. 
Fig. 2(a) shows the grid distribution on the symmetric plane. 
Because of the anticipated relatively steep gradients in velocity, 
pressure and temperature in zones such as at the chimney outlet, 
near the chimney wall and inside the chimney and collector, the 


Chimney 


a) grid distribution on the symmetric plane and (b) local grid distribution near the chimney. 
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grids need to be more concentrated than in the other zones, 
which is shown in the form as the quite darker areas. 
Fig. 2(b) displays the local grid distribution within the chimney 
and its outside ambience, where finer grids for the boundary layer 
near the chimney wall, collector wall and ground wall and the 
structured quadrilateral grids were adopted in order to reduce the 
grid number to improve the following computing speed without 
impairing the meshing quality. 


2.5. Computational procedure 


The computations have been performed by using the general 
purpose CFD program Fluent 6.3.26. The QUICK scheme was used 
to discretize the convective terms and a second order accurate 
treatment was used for the diffusion terms. The set of discretized 
algebraic equations was solved in a coupled manner. The itera- 
tions were continued until the relative error in the mass con- 
servation equation was below 5x107? and in the energy 
equation 1 x 1078. The simulation method using the relative static 
pressure in place of static pressure to analyze the whole pressure 
distribution of the whole system is the same as that being used by 
Pastohr et al. [89], Ming et al. [34,44,46] and Sangi et al. [40]. 

In order to test the grid-independent performance of the grid 
system selected in the numerical simulation, three test cases of 
the whole model under the same conditions (ambient crosswind 
is Om/s, and solar radiation is 857 W/m?) with grid numbers 
being 1,511,354; 1,674,272; 1,833,458 were tested. Numerical 
simulation results indicated that the volume flow rates of the 
chimney outlet are 789.19 m?/s, 822.77 m?/s, 833.35 m?/s, 
respectively corresponding to the three mesh systems listed 
above. By comparing the latter two mesh systems, we found that 
there was only a deviation of approximately 1.3% between these 
two results, which demonstrated the solutions in this case are 
grid-independent. The grid spacing and number of 1,674,272 is 
thus selected as the basic mesh system of this paper. 


2.6. Selection of ambient geometrical dimensions 


The juxtaposed two figures (Figs. 3 and 4) reflect the velocity 
contours in the simulated field. The difference lies in the size of 
the model, the right one is lengthened in both x and z direction 
relative to the left one. Due to the fact that there exists Karman 
vortex downstream, which would result in the periodic sway of 
the flow in the rear, it is plausible that the velocity contours 
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behaves unsteadily. Moreover, when the lower left quarter of the 
larger-scale model is cut off and analyzed, it is found that the 
contours is identical to those in the small one, thus testifying the 
reasonability of the results got from a model with just 200 m in 
the downstream domain. 

In order to find the impact of the scale of the simulated 
domain on the performance of SUPPS, based on the two different 
model, detailed analysis on the outflow updraft velocity under 
different ambient conditions, especially velocity, has been carried 
out. As can be seen from Fig. 4, a deviation lower than 3% can be 
detected from the changing curve of the two models shown in 
Fig. 3. Therefore, the power performance, which is the most 
important parameter we care about, can be considered as reliable 
results despite the smaller simulated scale. The appropriate main 
geometrical dimensions are thus shown in Fig. 1. 


3. Results and discussion 


For the SUPPS model whose geometrical dimensions mainly 
come from the prototype in Manzanares, Spain [1], the main 
factors that influence its characteristics are solar energy input 
into the system, ambient crosswind and fore-and-aft pressure 
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Fig. 3. Influence of ambient geometrical dimensions on system velocity distributions at G=857 W/m°. (a) L=200 m and (b) L=600 m. 
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drop across the turbine. In this paper, the solar radiation heating 
the ground surface of the collector can be regarded as heat flux, 
wind velocity magnitude at the top of the chimney and pressure 
drop across the turbine can be set beforehand. 

In this part, the effect of only the ambient crosswind on the 
performance of SUPPS was researched. A group of wind strength 
was adopted here: the ambient crosswind velocity at the height of 
the chimney outlet (200 m) changes from 0 m/s to 15 m/s, at an 
interval of 5 m/s. According to the logarithmic law of wind profile, 
in these conditions, the wind velocity at the height of 10 m is 0 m/s, 
2.2 m/s, 4.3 m/s and 6.5 m/s, tantamount to the Beaufort wind force 
scale 0, 2, 3 and 4, respectively. 

Generally, the solar radiation intensity in the deserts of north- 
west of China varies from 0 W/m? to 1400 W/m?. When the 
absorption of the soil and the energy loss being reflected as the 
form of radiation heat transfer back to the space are taken into 
consideration, the solar energy transferring to the air as thermal 
energy within the collector is just 60-70% of the total amount. 
Therefore, heat flux through the collector bottom ranges from 
0 W/m? to 980 W/m?, making 600 W/m? fairly representative of 
the common heating condition. Thereby, the heat flux through the 
ground of the collector was selected as being 600 W/m?, namely, 
the solar radiation is 857 W/m?, to analyze the effects of ambient 
crosswind on the velocity, pressure and temperature distributions 
of the system qualitatively. 
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3.1. Comparison of flow performances 


Fig. 5 displays the comparison of contours of velocity magni- 
tude at the symmetry plane in the whole simulated area when the 
velocity of the ambient crosswind at the height of 200m 
increases from 0 m/s to 15 m/s and the solar radiation intensity 
G remains 857 W/m’, constantly, where Uz99 m means the ambi- 
ent crosswind at the height of 200 m which is also the height of 
the chimney outlet. It is apparent that when no external wind is 
blowing horizontally, as shown in Fig. 5(a), the airflow in the solar 
chimney and the collector is axisymmetric: air nearby flows into 
the collector, accelerates gradually and runs centripetally to the 
foot of the chimney with energy being absorbed from the ground 
surface of the collector and air temperature increasing gradually; 
then, the heated air is sucked into the chimney and the updraft 
reaches its peak at a maximum speed of 19.3 m/s at the chimney 
bottom and finally passes through the chimney smoothly and 
steadily until it approaches the top outlet of the chimney. 
However, when the ambient crosswind at the top of the chimney 
increases to 5 m/s, as shown in Fig. 5(b), the negative influences 
are evident: the airflow at the chimney bottom is deflected with a 
maximum velocity magnitude less than 13.1 m/s; the outflow 
from the chimney outlet is also slanted downstream by the 
ambient crosswind instead of rising straightly upward. The same 
phenomenon occurs when ambient crosswind velocity at the 
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Fig. 5. Influence of ambient crosswind on velocity distributions in the symmetry plane at G=857 W/m?. (a) Uz99 m=0 m/s, (b) U200 m=5 m/s, (c) U200 m=10 m/s and 


(d) U200 m=15 m/s. 
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chimney outlet is 10 m/s shown in Fig. 5(c), namely, the peak air 
velocity inside the SUPPS decreases to 11.2 m/s with increasing 
ambient crosswind velocity. However, when approaching 15 m/s, 
as shown in Fig. 5(d), the ambient crosswind does not exacerbate 
the flow performance of SUPPS, instead, it seems to function 
oppositely: the peak velocity of airflow within the chimney 
bottom returns to 12.7 m/s, beyond our expectation according 
to the former downside influence of ambient crosswind. We can 
thus find that the ambient crosswind has a two-side effect on the 
SUPPS from the down-and-up performance of airflow velocity 
magnitude in the chimney. 

Fig. 6 contracts the static pressure when the ambient velocity 
differs from each other. From these figures, it is evident that the 
pressure undergoes a promotion though the chimney, no matter 
whether there is crosswind outside. Inasmuch no turbine is 
installed, wall friction mainly contributes to the depletion of 
pressure difference as the updraft flows inside. 

In order to further illustrate the phenomenon observed above, 
Fig. 7 shows the local velocity vectors near the chimney outlet at 
G=857 W/m? with different ambient crosswind velocities. 
To some extent, the vectors can reflect the velocity field as their 
length stands for the velocity magnitude and the arrow heads for 
the flow direction. As shown in these contrasted figures, the air 
from the chimney outlet flows straight upwards without ambient 
crosswind. Then, the outflow slants to the downstream with the 
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existence of ambient crosswind: the higher the ambient cross- 
wind velocity, the more deflective the outflow. It can be expected 
that, the weak ambient crosswind will inhibit the existent 
chimney outflow which is due to the buoyancy effect caused by 
the solar radiation on the collector ground surface, whereas the 
weak chimney outflow may be strengthened if ambient crosswind 
is strong enough. A very strong ambient crosswind horizontally 
flowing across the chimney outlet may produce a negative- 
pressure zone near the chimney outlet, resulting in an increase 
of the chimney outlet air velocity. 

Fig. 8 displays the local velocity vectors at the chimney bottom 
at G=857 W/m? with different ambient crosswind velocities. 
Comparing Fig. 8(a) to (b), we can easily see that airflow within 
the chimney abates with increasing ambient crosswind, whereas 
it rises with increasing ambient crosswind from the latter two 
figures. In addition, there are two phenomena to which we should 
pay close attention. On one hand, there is a vortex within the 
chimney bottom with the existence of the ambient crosswind, 
and it lies on the left side near the chimney wall as the ambient 
crosswind flows from the left to the right side. This vortex 
increases its scope with increasing ambient crosswind, leaving 
the updraft near the chimney wall on the right side being 
stronger. On the other hand, looking carefully at the airflow 
direction within the collector, we can find that the airflow sucked 
from the area, in the ambient crosswind’s downstream half of the 
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Fig. 6. Influence of ambient crosswind on pressure distributions in the symmetry plane at G=857 W/m2. (a) U200 m=0 m/s, (b) U200 m=5 m/s, (c) U200 m=10m/s and 
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Fig. 7. Influence of ambient crosswind on local velocity vectors near the chimney outlet at G=857 W/m?. (a) U200 m=0 m/s, (b) U200 m=5 m/s, (c) U200 m= 10 m/s and 


(d) U200 m=15 m/s. 


collector which in the figures is the half on the right, is weakened 
by the existing ambient crosswind. Especially, when the bottom 
ambient crosswind is strong enough, air may be blown to the 
outside of the collector instead of converging into the chimney. 
Therefore the enthalpy of the heated air is lost, which might be 
the main reason why ambient crosswind at low speed has 
negative effects on the performance of the SUPPS. 

Pretorius [90], in an attempt to find a “thermo-economically 
optimized” dimensional configuration of SUPPS by comparing 
construction costs to annual power output, tested many plant 
dimensions with a multitude of variables for different cost structures 
for SUPPS from 500 m to 1500 m high. His numerical model predicted 
plant configurations for which inflow of cold air into the top of the 
tower (called cold inflow) may exist in the chimney of the SUPPS, 
especially when air velocities through the plant are low (very slow 
upward air velocities). The probability of a plant experiencing cold 
inflow increases with decreasing collector diameter, chimney height 
and collector inlet height as well as increasing chimney diameter. 

In general, for Pretorius model, if the ratio H/D of the chimney 
height/chimney diameter is >5 there was no cold inflow, in the 
absence of crosswind, or with outside wind velocities of 2 m/s at 
10 m high and 4 m/s at the top of the chimney. In our model, with 
a smaller tower (Manzanares model), and higher crosswind 
velocities from 10 to 20 m/s inflow of air into the top of the 


tower seem to be experienced, as seen for the local velocity 
vectors in Fig. 8 for the bottom of the tower. 

Serag-Eldin analysis [84] revealed a huge degradation of 
performance with 10 m/s winds, and even with 2 m/s weak winds 
a considerable degradation occurs unless the collector inlet height 
is lower (7.5 m instead of 15 m). As a matter of fact, for his study 
he took the dimensions suggested by Haaf et al. [14,15] for a 
5 MW plant, but the collector entrance is too high (15 m). In our 
study this height is 2 m like in the Manzanares pilot plant. 

In these conditions, Serag-Eldin observed that the effect of the 
atmospheric crosswind is to blow the heated air downstream the 
collector, rather than up the chimneystack, thus reducing the air 
temperature in the chimneystack, the temperature difference 
between hot air in stack and atmospheric air outside is lower, 
the velocity of air in the stack is greatly reduced, and so is the 
flow rate. Moreover, since this motion is responsible for driving 
the wind turbine, the plant performance is severely impacted. 

In another work Serag-Eldin [91] addressed the problem of 
performance degradation of SUPPS when they are exposed to 
strong external atmospheric flows. In an attempt to reduce the 
hot air escaping sideways he proposed the introduction of con- 
trollable flaps (a quarter or a half circle perimeter) at the down- 
stream end of the collector, in order to reduce the proportion of 
hot air by-passing the chimney stack. 
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Fig. 8. Influence of ambient crosswind on local velocity vectors at the chimney bottom at G=857 W/m”. 


(d) U200 m=15 m/s. 


3.2. Comparison of relative static pressure contours 


Fig. 9 shows the relative static pressure contours near the 
chimney bottom at G=857 W/m?. Because this is the very place 
where airflow from the collector converges, the pressure under- 
goes the steepest gradient at the entrance of the chimney. Also, 
the negative relative static pressure, which represents the pres- 
sure difference between the airflow within the chimney and the 
stable atmosphere outside, has its minimum value at the turning 
from the collector to the chimney and then increases gradually 
through the chimney. By analyzing the contrast of pressure 
distributions as shown in Fig. 9(a-d), we can find that the 
minimum relative static pressures at the chimney bottom are 
— 254.37 Pa, — 118.36 Pa, —79.69Pa, —95.61 Pa respectively cor- 
responding to the four ambient crosswind velocities: Om/s, 5 m/s, 
10 m/s and 15 m/s. As the minimum relative static pressure of the 
SUPPS is a reflection of the system driving force analyzed by Ming 
et al. [34], the driving force of the SUPPS will also slump at first 
and then ascends with the increasing ambient crosswind velocity. 
Besides, the reversed flow field of model exposed in strong 
ambient crosswind is more expansive than in weak ambient 
crosswind; in particular, when external crosswind is at a speed 
of 15 m/s at the top of chimney, the vortex-shedding phenom- 
enon is pretty clear. Also, as shown in Fig. 9(c) and (d), the relative 
static pressure is positive in the downstream half of the collector 


(a) U200 m=0 m/s, (b) U200 m=5 m/s, (c) U200 m=10 m/s and 


when the ambient wind is relatively strong, which means that the 
pressure at this place is higher than that of the ambience at the 
same height. Thus the air within this zone where the relative 
static pressure is positive will flow from the collector inlet to the 
ambience which has been verified previously from the results 
shown in Fig. 8(c) and (d). 


3.3. Comparison of temperature contours 


Fig. 10 denotes the contrast of temperature contours in the 
symmetry plane of SUPPS exposed in different ambient crosswind 
velocities ranging from 0 m/s to 15 m/s. It is evident that the 
plume of the air outflow from the chimney outlet deflects more 
and that its average temperature decreases as the external cross- 
wind becomes stronger: the value of the average temperature is 
317.75 K, 303.84 K, 299.55 K, and 297.96 K for the four conditions. 
In addition, the plume scope of the air outflow influenced by the 
ambient crosswind becomes smaller and smaller with increasing 
ambient crosswind velocity. It is because on one hand, the 
temperature increase of the updraft decreases with increasing 
ambient crosswind velocity; on the other hand, convection heat 
transfer between the plume of the updraft and the ambience has 
been significantly enhanced with increasing ambient crosswind 
velocity. From what the comparison displayed, it is safe to draw 
the conclusion that ambient crosswind has a negative influence 
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Fig. 9. Influence of ambient crosswind on local pressure distributions at the chimney bottom at G=857 W/m’. 


and (d) U200 m=15 m/s. 


on the heat transfer process within the SUPPS, thus reducing the 
buoyant force originating from the density difference between the 
SUPPS and the ambience. It is apparent that the difference of 
temperature in the solar chimney among these conditions partly 
results from the imposed flow by the ambient crosswind under- 
neath the collector, as shown in Fig. 10. The crosswind heading 
towards the collector outlet in the downstream half of the 
collector brings about the outflow of heated air. This effect may 
contribute to the waste of buoyant force generated by heated air 
and deteriorate the SUPP performance. Furthermore, the external 
crosswind skating over the chimney wall and collector canopy 
surface may accelerate the heat transfer between the air within 
the SUPPS and the ambience, leading to the decrease of the 
updraft temperature increase, but this effect is comparatively 
insignificant. 


3.4. Comparison of system temperature increase, driving force and 
updraft velocity 


In this part, the numerical simulation results of the geome- 
trical model including SUPPS and ambience in Fig. 1 being given 
different solar intensity and ambient crosswind conditions are 
quantitatively compared and analyzed as shown in Figs. 11-13. 
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a) U200 m=0 m/s, (b) U200 m=5 m/s, (c) U200 m=10 m/s 


Fig. 11 displays the influences of external crosswind on the 
temperature increase of the updraft from the chimney outlet. As 
seen in this figure, the increase of updraft temperature decreases 
significantly with the increase of ambient crosswind velocity as 
long as the solar radiation intensity is over zero, and it increases 
with increasing solar radiation. Similarly, it is the ambient 
nethermost crosswind flowing into the collector that gives rise 
to the decrease of the updraft temperature increase. From the 
more integral data shown here, it is easy to observe that the 
SUPPS with higher solar radiation is more sensitive to ambient 
crosswind: the outlet temperature increase plummets from 
30.35 K to 13.51 K, about 17K drop, when solar radiation is 
1143 W/m? with ambient crosswind velocity at the chimney 
outlet increasing from 0 m/s to 5 m/s, compared with just 7 K 
drop when the solar radiation is 286 W/m2. This phenomenon 
might be mainly caused by increased cooling of the canopy and by 
the escape of air with relatively high enthalpy from the collector. 

Fig. 12 shows the influence of the ambient crosswind on the 
driving force of SUPPS with different solar radiations, where the 
minimum value of relative static pressure within the chimney 
bottom can directly reflect the driving force of the SUPPS. 

With solar radiation being more than 286 W/m? it looks like 
that with wind speed at 5 and 10 m/s the wind acts as a cover (cap 
effect) at the tower outlet which decreases the tower air flow. 
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Fig. 10. Influence of ambient crosswind on temperature distributions in the symmetry plane at G=857 W/m”. (a) U200 m=0 m/s, (b) U200 m=5 m/s, (a) U200 m=10 m/s 
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Fig. 11. Influence of ambient wind on the temperature increase of updraft. 


As shown in this figure, a very interesting phenomenon can be 
found as shown in this figure that, with solar radiation being 0 W/m?, 
the driving force of the SUPPS increases gradually with increasing 
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Fig. 12. Influence of ambient wind on the driving force of SUPPS. 


ambient crosswind velocity whereas it decreases significantly 
to valley values and then gradually increases with increasing 
ambient crosswind velocity when the solar radiation is over 
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Fig. 13. Influence of ambient crosswind on average velocity of updraft. 


0 W/m°. Apparently the ambient crosswind has positive influence on 
the driving force of SUPPS with solar radiation being 0 W/m?, the 
reason for this might be that lower velocity of ambient crosswind that 
directly enters the collector; higher velocity of ambient crosswind 
roars past the chimney outlet which will result in a negative pressure 
zone near the chimney outlet causing a suction effect of the SUPPS. 

The Bernoulli Principle might explain why the higher the 
velocity of the wind across the top of the chimney of the SUPPS, 
the faster it will draw air up the chimney (which in turn can 
increase power output). When wind blows perpendicular to the 
end of a tube, it creates a vacuum. The higher the wind velocity, 
the lower the pressure will be, and the more air will be drawn 
through an open-ended tube. As wind passes over the opening at 
the top of a chimney, the lower pressure sucks the air up the 
chimney. In other words with wind blowing with sufficient speed 
over the top of the chimney, the air pressure at the top is reduced, 
and tends to draw higher-pressure air from down below. Also the 
higher the wind speed, the higher the Venturi suction effect at the 
top of the chimney. The Bernouilli equations are simplifications of 
the Navier--Stockes Eq. (3). 

From this figure, when the ambient crosswind velocity is less 
than 10 m/s, it has negative influence on the driving force of 
SUPPS; and it will have positive influence on the driving force of 
SUPPS if it is higher than 10 m/s when the solar radiation is over 
0 W/m. This is because both negative and positive aspects of 
ambient crosswind influence the system simultaneously. Appar- 
ently, when ambient crosswind is relatively weak, the maximum 
pressure difference under these solar conditions differs much 
more from what they do when the ambient crosswind velocity 
approaches 25 m/s. This phenomenon indicates that with the 
ambient crosswind getting stronger, solar radiation functions less 
and is less crucial in determining the performance of SUPPS while 
the ambient crosswind becomes more and more predominant. 

Fig. 13 displays the relationship between the average velocity 
of the updraft and the ambient crosswind velocity under different 
solar radiation intensities. As shown in this figure, the air flow 
inside the chimney is notably accelerated with increasing ambient 
crosswind with solar radiation being 0 W/m?, which is in agree- 
ment with the result shown in Fig. 12. Similarly, there also exists a 
valley value of updraft velocity in each line with increasing 
ambient crosswind when the solar radiation is over 0 W/m7. Take 
the curve of G=1143 W/m? for example, we can see that at first 
the ambient crosswind undermines the average velocity of updraft 
and further impairs the performance of the SUPPS evidently. 


However, once the valley point is weathered, the positive effect, 
namely the suction effect of the SUPPS from the chimney outlet, of 
external crosswind plays a more important role. 


3.5. Influence of crosswind with turbine pressure drop 


Fig. 14 displays the influence of ambient crosswind on average 
velocity of updraft with turbine running when the solar radiation 
is 857 W/m7?, where the basic method for the control of the 
turbine pressure has been presented by Ming et al. [31,34] and 
Xu et al. [45]. As shown in this figure, when the turbine pressure 
drop remains constant, the relation between average velocity of 
updraft and the ambient crosswind velocity seems to be similar to 
that shown in Fig. 13. That is, the performance of the solar 
chimney experiences an overall down-and-up process and the 
minimum of average velocity of updraft occurs when ambient 
crosswind velocity is 10 m/s. Numerical simulation results with 
different solar radiations indicate that the similar phenomenon 
shown above can be found regardless of the pressure drop, at 
least in the scope of conditions the model of this paper is involved 
in. Meanwhile, the average velocity of the updraft decreases 
notably with increasing turbine pressure drop being given the 
same ambient crosswind velocity. This phenomenon further 
denotes that the ambient crosswind and the turbine pressure 
drop are two key factors that influence the performances of SUPPS 
independently. 

In order to evaluate the output power performance of SUPPS 
exposed to various ambient crosswind velocities, Fig. 15 displays 
how the output power of the turbine, which represents the ability 
of the SUPPS to convert the solar and wind energy into electric 
power, can be influenced by the pressure drop across it at 
G=857 W/m?. As for the curve computed in the absence of 
ambient crosswind, it is rather easy to see that it is so smooth 
that a parabolic curve of its overall profile can be deduced from the 
first half of it, similar to what has been observed by Xu et al. [45]. 
For this phenomenon, a convincing reason is proposed: according to 
Eq. (13), turbine output is determined by both turbine pressure drop 
and volume flow rate, only when the extent to which pressure drop 
increases is considerable can the rate of airflow decrease signifi- 
cantly, making the product of them to climb up at first and then run 
down with the pressure drop. 

In other conditions with ambient crosswind velocity being 
higher than 0 m/s, although the output power curves experience 
somewhat similarly, they yet differ much in specifics. Differences 


T H T T T T T T T T T 
10 + 4 
— —=— Ap Pa 
an a 
E —e— Ap=20 Pa 
a 8f —4— Ap =40 Pa | ] 
i 
E —— Ap,=60 Pa | 1 
S 6f 7 
= ] 
9 
2 
S 4p 5 
vo 
On 
S ] 
z 
2 2} : 
0 i i i i i 4 


0 5 10 15 20 25 
Ambient crosswind velocity (m/s) 


Fig. 14. Influence of ambient crosswind on average velocity of updraft with 
turbine running at G=857 W/m?. 


T. Ming et al. / Renewable and Sustainable Energy Reviews 16 (2012) 5567-5583 5581 


—#— Uj 99n=0 m/s 


40 + 


—2— Uj 99=5 m/s 
—4— Uj 99m=10 m/s 1 


—v¥—U 99 =15 m/s 4 


—+—U99m=20 m/s 


20 


Output power (kW) 
w 
© 


10 


0 40 80 120 160 200 
Pressure drop across the turbine (pa) 


Fig. 15. Influences of pressure drop on turbine output power. 


are the maximum of output power and the corresponding turbine 
pressure drop, as well as the protruding extent of the curves. 
When the ambient crosswind velocity is lower than 15 m/s, it has 
significant influence on the output power: the maximum output 
power can be over 50 kW at the turbine pressure drop of 160 Pa 
without ambient crosswind velocity whereas it is lower than 
10 kW with ambient crosswind velocity being 10 m/s. It can also 
be seen from the three output power curves (U200 m=0 m/s, 5 m/s, 
10 m/s) that their changing trend is very similar. Thus the ambient 
crosswind has negative influence on the system output power. 
However, the output power curves of U299 m= 15 m/s and 20 m/s 
as shown in Fig. 15 differ greatly from the curves with ambient 
crosswind velocity being lower than 15 m/s. Take the curves of 
U200 m=20 m/s, 15 m/s and 10 m/s as an example, we found that 
being given certain turbine pressure drop, the system output power 
increases with increasing ambient crosswind velocity, which is 
notably different from the changing trend shown with ambient 
crosswind velocity is lower than 15 m/s. Thereby, the ambient 
crosswind has positive influence on the output power of SUPPS. In 
addition, the maximum output power of SUPPS can be reached at 
certain turbine pressure drop point with ambient crosswind velocity 
being lower than 15 m/s, whereas it can be reached during a broad 
scope of turbine pressure drop with ambient crosswind velocity 
being higher than 15 m/s. For instance, when the ambient crosswind 
velocity is 20 m/s, the system output power is about 13-14 kW with 
the turbine pressure drop changing from 60 to 160 Pa. 


4. Conclusions 


Through the detailed numerical simulations, this paper has 
investigated the influences of various ambient crosswind velo- 
cities on the performance of SUPPS. The solar radiation and 
turbine pressure drop were all being taken into consideration in 
the analysis. The geometrical model including the SUPPS of the 
Spanish prototype and its ambience were developed together 
with the mathematical models describing the fluid flow, heat 
transfer and output power performances of the system. The fluid 
flow, temperature and pressure distributions as well as the 
temperature increase, driving force and output power parameters 
were also presented and analyzed. 

The numerical simulation results revealed that the ambient 
crosswind has significant influence on the fluid flow, heat transfer 
and output power performances of the SUPPS. The feature of the 
influence has two sides, both positive and negative. 


The results have indicated that at low speed the wind effects 
can be attributed to 


e the thermal flow effect at the level of the greenhouse canopy 
(wind increases convective heat losses from the collector roof 
to the environment), 

e the entering air at the bottom of the tower causing air flow 
distortions and resistances, the reduction of pressure or 
temperature difference between the inside and outside of the 
tower, 

e non-uniform air temperature distribution inside the tower and 
flow distortions, 

e the wind flow covers the exit and deflects the plume. As a 
result the reduction of the effective area of the chimney outlet 
causes a throttling or a cap effect which decreases the output 
air flow, 

e possible infiltration of wind into the tower through its top 
resulting in choking flows inside the tower. 


The results seem to indicate that at high speed the wind effects 
can be attributed to 


e a wind suction effect (Bernoulli type, might be increased by a 
Venturi deflector) at the tower outlet which increases the 
tower airflow. 


It was found that the SUPPS’s performance would be deterio- 
rated if the ambient crosswind velocity was below 15 m/s, and 
the performance was gradually improved with ambient cross- 
wind velocity being higher than 15 m/s. The trend under different 
external ambient crosswind conditions may contribute to the 
optimal design and control for the running of turbine of SUPPS. 
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